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Can we measure the accretion efficiency of Active Galactic Nuclei? 
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ABSTRACT 

The accretion efficiency for individual black holes is very difficult to determine accurately. 
There are many factors that can influence each step of the calculation, such as the dust and 
host galaxy contribution to the observed luminosity, the black hole mass and more importantly, 
the uncertainties on the bolometric luminosity measurement. Ideally, we would measure the 
AGN emission at every wavelength, remove the host galaxy and dust, reconstruct the AGN 
spectral energy distribution and integrate to determine the intrinsic emission and the accretion 
rate. In reality, this is not possible due to observational limitations and our own galaxy line 
of sight obscuration. We have then to infer the bolometric luminosity from spectral measure- 
ments made in discontinuous wavebands and at different epochs. In this paper we tackle this 
issue by exploring different methods to determine the bolometric luminosity. We first explore 
the trend of accretion efficiency with black hole mass (e oc M""'^-^) found in recent work 
by Davis & Laor and discuss why this is most likely an artefact of the parameter space cov- 
ered by their PG quasar sample. We then target small samples of AGN at different redshifts, 
luminosities and black hole masses to investigate the possible methods to calculate the accre- 
tion efficiency. For these sources we are able to determine the mass accretion rate and, with 
some assumptions, the accretion efficiency distributions. Even though we select the sources 
for which we are able to determine the parameters more accurately, there are still factors af- 
fecting the measurements that are hard to constrain. We suggest methods to overcome these 
problems based on contemporaneous multi-wavelength data measurements and specifically 
targeted observations for AGN in different black hole mass ranges. 
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1 INTRODUCTION 

Accretion onto a supermassive black hole is the current accepted 
mechanis m to expl ain Active Galactic Nuclei (AGN) activity 
(iSalpeteril964 iLvnde n-Bell 1969). This extreme process is respon- 
sible for radiative emission on a wide wavelength range, with most 
of the intrinsic power emitted in the optical, ultra-violet (UV) and 
X-ray wavebands. With multi-wavelength data it is possible to de- 
termine the AGN spectral energy distribution (SED), and investi- 
gate the physical mechanisms behind the radiative emission at each 
energy. Earl y work set the standard model for a thin black hole ac - 
cretion disc ('Shak ura & Sunvaev|[T973l : lNovikov & T home" 1 973'). 
and s howed that it could explain the observed emission (Shields 
Il978t) . establishing a relation between the SED peak in the UV re- 
gion and the thermal emission from the accretion disc. This pic- 
ture has been explored with the development of more complex 
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mod els, and although there are some observational discrepancies 
(see lKoratkar&B"lae3ll999 ' for a review), the thin accretion disc 
paradigm remains the standard adopted model. 

The accretion parameters such as the rate of accreted mass 
(M) and efficiency in converting mass to radiative energy (e = 
L/{M(?)) are closely related to the energy output and play a deci- 
sive role on AGN behaviour and evolution. The study of their vari- 
ation with time, environment and type of source can give insights 
into black hole accretion physics and constrain the current models 
for the evolution of AGN properties. Soltan ( 1982) showed that it is 
possible to constrain the average efficiency in converting accreted 
mass to emitted radiative energy (accretion efficiency hereafter), 
by comparing the integrated AGN luminosity with the total mass 
accreted onto black holes. Later work used a related method to de- 
termine the average accretion efficiency by measuring the X-ray 
background jFabian & Iwasawalll999l : lElvis et al. I l2002l) . In stud- 
ies of AGN evolution, where the redshift dependent AGN luminos- 
ity function is compared with the local black hole mass function 
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2OO9I). the accretion efficiency is one of the free parameters as- 
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sociated with the accretion process. These can be determined on 
average for AGN, but when taking the observed values for Ed- 
dington ratio (Lboi/^Edd) and obscu ration into account, the pos - 
sible parameter range is very wide JRaimundo & FabiMl l2009h . 
More observational results are needed to constrain the evolution 
of black holes and AGN properties. Direct measurements of indi- 
vidual accretion parameters would allow to constrain and modify 
the models to better take into account the variety of AGN stages 
and environments. Recent work explored the possibility of deter- 
mining the mass accretion rate and accretion efficiency for indi- 
vidual AGN, using accretion disc t heory argum ents and available 
data on the emitted radiation. . Collin et alj ( l2002h used a simple ac- 
cretion disc model to determine the mass accretion rate from the 
observed optical monochromatic luminosity, and explored in de- 
tail the uncer t aintie s associated with the simple disc assumption. 
iBian & Zhaj ( I2OO2I) . adopted a similar relation to determine mass 
accretion rates and accretion disc inclinations, and used this as a 
starting point to derive accretion pr operties for radi o loud and ra- 
dio qu iet AGN jBian & Zhaoll2003h . More recently. [Pavis & LaoJ 
( I2OI ll) (DLl 1 hereafter), compared the simple accretion disc model 
with more sophisticated ones and determined individual mass ac- 
cretion rates and accretion efficiencies for 80 PG quasars. Their ac- 
cretion efficiency values showed an increase with increasing black 
hole mass, and a wide distribution along the allowed range for the- 
oretical black hole spin predictions. From accretion disc theory, 
the accretion efficiency is related with the spin of the black hole 
and the procedure described above provides a method to measure 
individual black hole spins. The relation between the two is set 
in first order by accretion disc theory, but note that effects such 
as magnetic torques can increase the amount of energy released 
in the disc and the relation betwe en the spin and accretion effi- 
cienc y becomes more complex ("e.g. lGammi ell999l : lAgol&Krolikl 
I2OOCI) . The spin distribution can provide clues on the history of 
supermassive black holes, with low values of spin associated with 
black hole growth by random minor mergers and high values as- 



socia t ed with growth by continuous gas accretion 1 Volonteri et al.l 



I2OO5I : iKing & Pringlel2006l : iBerti & Volonterill2008l) . The spin dis- 
tribution determined from this method can also be compared with 
alternative approaches to measure black h ole spin, such as spectral 
analysis of the X-ray Fe Kg line emission 1 Fabian et al .ll989l:lLaoJ 
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Figure 1. Example of two different spins for the accretion disc model of 
IRXS J045205.0+493248. The lower pair of lines (blue) represent an ac- 
cretion disc with spin of a = 0.1, while the upper lines (black) represent 
a maximally rotating black hole. Solid and dash-dotted lines are the SED 
before absorption for each spin case, the dashed and dotted lines are the ob- 
served SED (after absorption). The vertical line shows the B-band observed 
wavelength, which is relatively insensitive to changes in spin. 



2 MODEL 

In this section we present an overview of the model used to deter- 
mine the mass accretion rate and accretion efficiency in this work. 
The method makes use of the fact that the AGN spectral energy 
distribution at optical wavelengths is relatively insensitive to black 
hole spin, which allow us to establish a relation between the spe- 
cific luminosity and mass accretion rate (see Fig.[T]but also Fig. 1 
of DLll). The mass accretion rates are determined from assuming 
a simple accretion disc model, the measured flux at 4392 A and 
an estimate for the black hole mass and disc inclination. The flux 
emitted by the accretion disc is the result of the integrated black- 
body emission (B^) at each radius R: 



2n cos(j) 



B^[T^fi)RdR 



(1) 



Assuming isotropic emission and an effective temperature de- 
pending on the black hole mass, mass accretion rate and radius 
(Tefi oc 7? 4 ), the relation between the luminosity at an optical 
frequency (L^) and the mass accretion rate (Af ) is: 



In this work, we discuss the uncertainties involved in the de- 
termination of accretion efficiencies and investigate the correlation 
between black hole mass and accretion efficiency. In Section 2, 
we present the main accretion disc model and how the accretion 
properties are determined. In Section 3, we analyse the correla- 
tion found by DLll and describe how the sample selection effects 
decrease the allowed parameter space and produce the correlation 
found in their work. In the last sections we present several methods 
to determine the bolometric luminosity and the accretion efficiency 
of AGN based on observational data. We analyse two samples of 
low and higher redshift AGN and discuss our results in Section 5. 
Section 6 summarises our main conclusions. Throughout this pa- 
per we adopt standard cosmological parameters of Ho = 71 km 
s"^ Mpc"\ ri,n = 0.27 and Q.^ = 0.73. 
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Where M is the black hole mass, i the disc inclination and x — 
hu /(kBTe f f). This relation is t he same as the corrected version 
of DL 1 1 in lLaor& David ( l201lh . The integral can be solved with 
Gamma and Zeta functions which give a numerical value of 1.93. 
Changing units, we obtain: 



M = 1.53 



10-*5 cos(i) 



3/2 



10** 



(MQ/yr) 



(3) 



with M in Mq and vL^, (or Lopt) in erg s . This equation is de- 
fined for the wavelength used in this work (4392 A), but the de- 



pendence in frequency is M oc {lyLi, 



,3/2 



. The accretion effi- 
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Figure 2. Black hole mass and optical luminosity distribution for the 
sources in the DLll PG quasar sample (red crosses). The black dashed 
lines represent the limits for our simulated sample. They are defined by flux 
and volume limitations and Eddington ratio biases. The dashed vertical grey 
lines mark the redshift range probed by the sample. 



ciency can then be determined from the mass accretion rate: 



Figure 3. Simulated distribution of accretion efficiency and black hole mass 
for AGN (grey circles) based on the flux and volume limits from Fig.|2]and 
a constant bolometric correction of 9.1. Top panel: Comparison with the 
DLll observed sample (red crosses). Bottom panel: Lines indicate qualita- 
tive behaviour or our sample's distribution. The green solid lines indicate 
the effect of luminosity limits. Blue dashed lines indicate the general effect 
of two Eddington ratio limits. The parameter space possible to probe with 
these conditions is indicated by the distribution of grey circles. 



Mc2 



(4) 



but this imphes having an accurate measurement of the observed 
bolometric luminosity (Lboi)- Here we decide to use the isotropic 
bolometric luminosity Lboi = A-kD^ x F, where F is the flux. For 
a disc, the bolometric luminosity should be a factor of 2 x cos(i) 
lower for a simple case and have a few extra factors for the rela- 
tivistic case. We nevertheless keep this formulation throughout the 
paper for a clear comparison between different methods of deter- 
mining Lboi- 



3 CORRELATIONS AND SELECTION EFFECTS 

The approach described in the previous section was used by DLl 1 
to determine the mass accretion rate and accretion efficiencies for 
a sample of 80 PG quasars. They explored complex models of ac- 
cretion discs to determine the best relation between the luminosity 
measured at an optical wavelength (4861 A in their case) and the 
mass accretion rate. The bolometric luminosity for each of these 
quasars was determined by integrating the spectral energy distribu- 
tion inferred from available data, interpolating between measured 
values and extrapolating to the regions where measurements are 
non-existent. In this section we use the DLl 1 best-fit relations and 
explore the correlation they found between the accretion efficiency 
and the black hole mass. DLl 1 discuss the selection effects in their 
sample, here we analyse these effects in more detail and simu- 
late the resulting allowed parameter range. From the discussion in 
the following paragraphs, it will become clear that the correlation 
arises from the source selection criteria. After this analysis we sim- 
ulate a test sample with uniform accretion efficiency and study the 
effect of uncertainties in inclination and black hole mass measure- 
ments. 



3.1 Sample selection 

The 80 PG quasars do not constitute a complete sample. Initially, 
these quasars were selected as part of the Palomar Bright Quasar 
Survey (BQS), a sub-sample of the Palomar-Green Survey (PG) of 
stellar-like objects with ultraviolet excess w hich is a volume and 
flux limited survey dSchmidt & Greenlll983h . The selection crite- 
ria included the presence of broad lines, a significant redshift (with 
the closest quasar at z = 0.025), and most of the emission con- 
centrated in the inner regions of the source (star-like appearance). 
DLll used the spectrosco pic data on a red s hift l imited sample 
(z < 0.5) of the BQS from boroson & GreerJ ( Il992h . In summary, 
the sample is biased to luminous quasars with high Eddington ra- 
tios, and subject to flux and volume limitations. The low Eddington 
ratio sources are likely to be missed from this sample, due to the 
selection criteria based on the star-like appearance as mentioned 
in DLl 1. For the low Eddington sources, the host galaxy luminos- 
ity can dominate the emission in the optical bands and most of the 
radiation will not come from the central source, failing the star- 
like selection criterion. The sample described above was neverthe- 
less very suitable for the study: the multi-wavelength data available 
makes it easier to determine the bolometric luminosity and calcu- 
late the relevant parameters for a large sample. This approach is 
valid for individual objects, but one needs to consider carefully the 
biases when determining general trends and correlations. 

From the analysis in DLll, the trends found using the 
more complex accretion disc models for a constant inclination of 
cos (i) — 0.8 were: 



M = 3.5 Mq yr~ 
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Figure 4. Simulated sample in the same conditions as in Fig. [5] with over- 
plotted line limits to indicate luminosity and Eddington ratio limits. Three 
panels illustrate the effects of different changes. Top panel: Effect of includ- 
ing lower Eddington ratios, lower limit in this panel is A > 0.001. Middle 
panel: Include sources at higher redshift with higher luminosity. Bottom 
panel: Lower the flux limit for source detection. 
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Figure 5. Simulated sample (grey circles) with DLll PG quasar sample 
over-plotted (red crosses) to illustrate the effect of using different bolomet- 
ric corrections. Top panel: higher bolometric correction (Xopt = 20). Bot- 
tom panel: Simulated distrib ution when con sidering a bolometric coirection 
weakly dependent on Lopt jMarconi et al.li2004) . In this case the values of 
accretion efficiency show a larger spread. 



Assuming a bolometric correction Kopt = iboi/-t'opt, the accre- 
tion efficiency will be: 

e = 0.063 7.-f (7) 

In DLl 1 the PG quasars are distributed along a region in the e vs 
Mbh that corresponds to a trend of e cx M^'^. Here we explore this 
correlation keeping in mind the characteristics of the sample and its 
selection criteria. In Fig.|2]we show the distribution of optical lumi- 
nosities (-Lopt) and black hole masses {AIbh) as a function of red- 
shift for the DLll sample. The redshift range is determined from 
the constraints discussed above, being limited to z ~ (0.025 - 0.5). 
We argue that these constraints, together with an almost constant 
bolometric correction, shape the distribution of accretion efficien- 
cies for these sources. The correlation found between e and A/bh in 
DLll is caused by the narrow range of parameters that are probed 
by this sample and not by a specific dependence of accretion effi- 
ciency with mass. If we take the distributions in Fig.|2]and use Eq.|7] 
we can only probe the parameter range indicated by the grey circles 
in Fig. [3] and therefore obtain a trend between efficiency and back 
hole mass. To verify this argument we generate a uniform sample 
of AGN with properties similar to the ones in DLll. First we start 
by determining the Lopt distribution for a flux limited sample with 
uniform distribution within the limits given by the dashed lines in 
the top panel of Fig.|2] black lines for the luminosity limits and grey 
lines for the redshift range. The limits we used can be understood 
qualitatively: at each redshift, the range of optical luminosities is 
determined on the lower limit side by the minimum flux for selec- 
tion, and on the upper limit by both the volume limit of the survey 
and the shape of the luminosity function, which predicts a decrease 
in the num ber volume of sou rces with increasing luminosity (see 
for example fSovle et alj200Cll for the _B-band luminosity function). 
It is nevertheless difficult to determine the flux limit for this sam- 
ple, Schmidt & Greeri (1983) give a _B-band magnitude limit for 
each exposure field, but the values of Lopt are measured at 4861 
A, and we would have to assume a spectral shape to determine 
the frequency specific flux limit. Instead, we assume a minimum 



uFi, = 3 X 10~^'^erg cm~'^ s~^ for every case, which is lower than 
the average _B-band mag nitude of 16.16 and spectral index of -0.5 
jSchmidt & Greenlll983h . but reflects well the lower limit observed 
in DLll. For the upper limit, we did not consider the luminosity 
function distribution (or volume limitations), which would be be- 
yond the scope of this paper, but assume a flux limit type border that 
contains almost all the sources observed in DLll. For the bottom 
panel of Fig.|2l the mass distribution in DLll corresponds to the 
virial estimated black hole masses and shows a broader distribution. 
To simulate this sample we generated redshift constrained sources 
within a lower and upper limit range, also represented as black and 
grey lines in the figure. Qualitatively, these limits can be understood 
as the lower and higher possible black hole masses at each redshift 
if we consider the luminosity limits, a bolometric correction i^opt 
(see discussion below) and the typical Eddington ratio of the DLl 1 
sample (0.04 < A < 1.0), with A = Lboi/(1.3 x 10^**Mbh/Mq). 
Using the simulated sample within the limits of Fig. |2] we then cal- 
culate the accretion efficiency for every source using the best fit 
relations from DLll shown in Eqs. |5]to|7] In Fig. [3] we present 
our simulated sample in grey compared with the DLll results in 
the top panel, and explain qualitatively the boundary limits in the 
bottom panel. This simple sample shows the effect of redshift, lu- 
minosity and mass limits in limiting the parameter space probed. In 
the bottom panel we plot in solid green the constraints due to lower 
and upper luminosity limits and in dashed blue the absolute limits 
due to the Eddington ratio constraints. In Fig.|4]we illustrate what 
would happen if lower Eddington ratio objects are considered (up- 
per panel), higher redshift sources (and consequently higher Lopt) 
are included (middle panel) and if a lower flux limit is consid- 
ered (bottom panel). We over-plot some of the lines from Fig. |3] 
as a guide. DLll find an almost constant ratio of Lboi to Lopt, 
in agreement with the low dependence of this ratio with bolomet- 
ric luminosity ( Marconi et al . 2004). In our calculations we assume 
a constant ratio of 9.1 which is the average found by DLll. The 
panels described above show that even with a constant bolometric 
correction, the parameter range will constrain the type of correla- 
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Figure 6. Simulated distribution of accretion efficiency and black hole mass 
for a sample with a similar parameter range to DLll (grey circles). The 
parameter range constrains the values of accretion efficiency and the corre- 
lations one can find. Our simulated sample shows a trend of e oc Af"'^® 
similar to the e oc M"'^ found by DLll. For a different parameter space 
the correlation would be different. 

tion one would find between the accretion efficiency and the black 
hole mass. If we include a higher bolometric correction or a spread 
in values, the accessible parameter range would also be affected. In 
Fig. [5] we plot how our sample (grey circles) would be distributed 
in the plane compared with the DLl 1 sources (red crosses). In the 
top panel we assume a higher bolometric correction (A'opt = 20), 
and in the bottom panel, t he effect of considering the Lboi depen- 
dent iM^onTet^] l l2004h bolometric correction. The first condi- 
tion would cause a shift, and the second a larger spread in e. 

The procedure described above yields a trend of efficiency 
with black hole mass by simply constraining the Lopt, A/bh and 
redshift distribution to be similar to DLl I (Fig.[2}- In the next para- 
graph we add more specific constrains such as the Eddington ratio 
range and determine the slope determined from the simulated e vs 
A/bh. 

We simulate a sample to explore the correlations one would 
see simply based on the parameter space available to the DLll 
quasar sample. We use the constraints mentioned above for lumi- 
nosity and black hole mass as a function of redshift and use a con- 
stant bolometric correction. To better represent the typically high 
Eddington ratio quasars in DLll, we add a A > 0.04 cut-off but 
leave the possibility of A > LO to include some of the sources 
which, from their calculations, show an Eddington ratio above one. 
We also limit the mass range to the PG quasar minimum and maxi- 
mum values log(A/BH/MQ) = (6.0 - 9.6) and log (Lopt [erg s~^]) 
= (43.4 - 46.5). Using Eq.|7] we determine the accretion efficiency 
and plot our results in Fig. |6l We over-plot our simulated sample 
(grey dots) and the results by DLl 1 (red crosses). We can see that 
with a constrained distribution of Lopt and A/bh due to selection 
effects, we can only probe a small range of parameters. Our sim- 
ulated sources show a trend e oc M^'^^ similar to the one found 
in DLll (e oc Af°'^), illustrating that this is mainly caused by the 
allowed parameter space covered by the sample. In general, if our 
AGN has a lower Eddington ratio and/or a lower Lopt, it should 
populate the empty upper region of the plot. If we select sources 
from larger volume surveys, we should also be able to fill in the 
region of high mass and high luminosity accretion efficiencies. The 
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Figure 7. Intrinsic distribution of accretion efficiencies and black hole 
masses for our test sample. Grey region indicates area of efficiencies rang- 
ing from the minimum counter-rotating black hole (e = 0.035) to a maxi- 
mally spinning black hole (e = 0.31). Black dots represent generated uni- 
form sample of AGN selected to have Eddington ratio in the range (0.01 - 
1.0). 

trend between accretion efficiency and black hole mass would con- 
sequently be different for such a sample. This will be visible in later 
sections when we analyse both local and high redshift sources and 
locate them in this plane. 

In terms of the actual values of accretion parameters, some 
sources show very low e and others Eddington ratios above one. 
The simulated sample discussed above also shows these properties 
because we did not constrain it further than the DLl 1 distributions. 
As discussed in DLl I, there are factors not included in their calcu- 
lations that would raise the values of accretion efficiency obtained, 
namely the host galaxy contamination and the dust reddening. The 
host galaxy contribution increases the value of Lopt , correcting the 
measured values for this effect would move the red sources to lower 
Lopt, which means up and to the left in Fig.|6] The dust correction 
acts to decrease the observed bolometric luminosity. Correcting for 
this effect would alter the parameters in the same way that an in- 
crease in the bolometric correction would do, which means moving 
the sources to higher accretion efficiencies as can be seen in the ex- 
treme case illustrated in Fig.|5] upper panel. As mentioned in DLl I, 
it is unlikely that e is constant for all sources, simply based on the 
parameter limits plotted in our figures, we can see that at least the 
Eddington limit will cause some of the regions to be free of sources. 
The other trends will always depend on the redshift range probed, 
sensitivity limits and the luminosity and black hole mass functions. 



3.2 Uncertainties in inclination and black hole mass 

In this section we simulate another test sample to investigate the 
effects of uncertainties in the inclination and black hole mass in the 
accretion efficiency calculations^ We start by simulating a sim- 
ple set of AGN with well established intrinsic properties: uniform 
accretion efficiency in the theoretical range 0.035 < e < 0.31, 

^ Please refer to Section 4.2 of DLl 1 for their discussion on the uncertain- 
ties in A/bh ■ 
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Figure 8. Intrinsic optical luminosity (vLy ) as a function of black hole 
mass calculated for a sample with accretion efficiency plotted in black in 
Fig. I?] See Section l331 for calculation details. 



black hole mass 10'' < Mbh /Mo < 10^°, Eddington ratio 
0.01 < A < 1.0, constant bolometric correction A'boi = 9.1 and 
inclination values between 10° and 80° uniformly distributed in 
cos (i) to simulate the solid angle weighting. These properties con- 
strain the Lopt distribution for these sources. Our aim is then to 
infer if, starting from this intrinsic Lopt distribution and using Eqs. 
[3] to [7] we would be able to retrieve the real properties of the sam- 
ple. 

In Fig. |7] we show the accretion efficiency uniform distribu- 
tion before (grey region) and after applying an Eddington ratio cut 
(black dots). In Fig. [8] we plot the intrinsic L^pt distribution calcu- 
lated using Eq.|7]and the distribution of cos(i) (the implicit depen- 
dence in this equation is e oc (cosi)^^'^). We then do the backward 
calculation and use this Lopt distribution and Eq. |7] to determine 
the measured accretion efficiency. The difference this time is that 
we use an observed black hole mass constrained to be within a fac- 
tor of three of the intrinsic black hole mass (i\/intr/3 < M^hs < 
3 X Mintr) to simulate the typical observational uncertainties in 
measuring this property. The inclination angle in taken as a constant 
value of cos(i) = 0.8, although the intrinsic distribution includes 
angles between 10° and 80°. The inferred accretion efficiency dis- 
tribution is then calculated using Eq.|7]and the results are plotted 
in Fig.|9] The diagonal cutoffs depend on the Eddington range as- 
sumed, and the diagonal spread to higher and lower values of ac- 
cretion efficiency are caused by the uncertainties in black hole mass 
and the effect of assuming a constant angle to a sample that has an 
intrinsic spread in inclination. With this test we can see that values 
below or above the predicted accretion efficiency limits can in part 
be due to uncertainties in the black hole mass and inclination angle. 
This test also shows that a uniform accretion efficiency distribution, 
when selected as a function of Eddington ratio, will show a trend 
with black hole mass. Although the correlation found by DLl 1 is 
likely due to the parameter space probed, their work has clearly set 
the method to determine accretion efficiencies for individual AGN. 
The parameter space covered by their sample would necessarily 
cause a distribution of sources in a constrained e vs A/bh region. If 
we consider other samples and their physical properties, based on 
the DLl 1 simple arguments, we can now calculate the accretion ef- 
ficiency for AGN. As referred to in the previous subsections, there 
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Figure 9. Distribution of intrinsic (black dots) and inferred (green triangles) 
accretion efficiency and black hole masses for a test sample of AGN. The 
inferred values assume a factor of 3 error in the black hole mass and a fixed 
inclination. This causes a spread in the diagonal direction creating a differ- 
ence between what we infer from the data and the intrinsic distiibution. 



are still issues in measuring these properties which can affect our 
results. We mentioned the effect of host galaxy correction, inclina- 
tion and black hole mass, but there are more factors that can cause 
problems in these calculations. In the next sections we will discuss 
the method of determining mass accretion rates and accretion effi- 
ciencies, the possible uncertainties involved in the calculations and 
apply this to examples of AGN with different characteristics. 



4 ACCRETION EFFICIENCY MEASUREMENTS 

In the next sections, we will explore several methods to estimate 
the accretion efficiency of individual AGN and discuss the uncer- 
tainties associated. We will use the standard accretion disc relations 
described in Section|2] Although more complex disc models exist, 
such as the ones explored in DLll, we use this simple approach 
because it reflects the main physical dependences in the accretion 
parameters, while allowing to evaluate analytically the source of 
uncertainties. There are two main steps in the accretion efficiency 
determination: the calculation of the mass accretion rate and the 
determination of the bolometric luminosity. The former can be cal- 
culated from the emission at optical wavelengths, using Eq. |3]and 
an inclination value. If we assume that this equation reproduces the 
trend between optical luminosity and accretion rate, the main un- 
certainties are the disc inclination, the black hole mass, the host 
galaxy contribution to the emission, and the presence of dust. As 
can be seen in Fig. [T] AGN with similar values of mass accretion 
rate, can differ in their bolometric luminosities due to the effect of 
spin: a source with higher spin will have a higher bolometric lu- 
minosity. To constrain the accretion efficiency, and consequently 
the spin, one needs to have an accurate estimate of the bolometric 
luminosity. This has always been a complicated issue, AGN emit 
radiation over a wide spectral range due to different emission mech- 
anisms. Ideally, knowing the spectral energy distribution would al- 
low us to determine the bolometric luminosity, but in practice, this 
is constrained by the fact that for most sources, the wavelength cov- 
erage is not very good, and that there are spectral regions impossi- 
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Figure 10. Observed frame spectral energy distribution example of a mod- 
elled turnover for one of the SDSS sources. Black circles indicate the op- 
ticalAJV data points corrected for galactic reddening and the X-ray data 
point. Dashed line is the model for observed emission and solid line is the 
model for intrinsic emission (before galactic absoiption). 



ble to probe, such as the far-UV regions where emission is obscured 
by the gas and any intervening dust in our galaxy. 

Here we explore possibilities of determining the bolomet- 
ric luminosity for different AGN samples. One is, of course, the 
method adopted by DLll. They gathered all data points and spec- 
tral slopes available on the PG quasars and assumed typical quasar 
slopes for regions without data coverage. Their SEDs are based 
simply on the data and do not attempt to follow a standard AGN 
SED model. This allows freedom in the SED shape and variations 
for individual AGN, but as mentioned in their work, can also re- 
sult in extreme spectral shapes for a small sub-sample of sources. 
Another approach is to select objects that show a turnover in their 
opticalAJV data, assume that it constrains the spin (see Fig.llOt. and 
integrate the optical to X-ray part of the SED to obtain the intrinsic 
bolometric luminosity. In the current AGN paradigm, we expect to 
see reprocessed emission in the infrared (IR), due to the optical/UV 
disc emission being re-emitted by material within the surrounding 
parsecs of the AGN. We use this alternative argument to determine 
the bolometric luminosity from the integrated IR emission and the 
X-rays, keeping in mind that geometry also plays a role in this case, 
since not all of the radiation is reprocessed. 



4.1 Optical/UV and X-ray simultaneous observations 

Variability is one of the important factors to take into account when 
reconstructing an AGN SED. It is well known that most AGN show 
some sort of flux variability, with timescales that depend on the 
waveband observed. We first start by determining the mass accre- 
tion rate for a sample of objects with simultaneous observations in 
the optical/UV and X-rays. 

The sample consists of 26 local AGN (redshift up to ^ 0.08) 
from the 5w;/f/BAT 9-month catalogue with observations in the 
UV-optical (UVOT) and X-ray telescopes (XRT). These sources 
were selected due to their low absorptio n and low spectral com - 
plexity and have been analysed in detail in lVasudevan et al] j2009l) . 
on a study of their physical properties and bolometric corrections. 
The UVOT flux values we use are from the mentioned work as 
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Figure 11. Example of the SED of one of the excluded sources from our 
analysis. Observed data points show a shape that can indicate the presence 
of dust. The solid line shows the best-fit intrinsic SED for this sourc e when 
including a fixed value of E(B - V) = 0.25 jwinter et al]|201Q) . Not 
considering the contribution of dust can underestimate the bolometric lumi- 
nosity. 



well, and have been corrected for host g alaxy emission and Galac- 
tic extinction (see sections 3.1 and 3.2 of lVasudevan et al. 2009 for 
details on this procedure). For most sources there a re six data points 
corre sponding to the UVOT filters (see Table 1 of lVasudevan et all 
l2009h . which we use in this work to determine directly the mass 
accretion rate. Our final sample consists of 22 sources, since 4 of 
the 26 sources show clear signs of SEDs affected by dust and were 
therefore excluded from our analysis: MCG-06-30-15, Mrk 590, 
Mrk 766 and ESO 490-G026 (see Fig.[TT]as an example of dust in 
SED). 

The optical monochromatic luminosity {Lv) is measured at 
the rest-frame wavelength of 4392 A. For the AGN at low redshift, 
since the observed and rest-frames do not differ by much, we use 
the B-band flux value at 4392 A observed wavelength and extrapo- 
late it linearly to obtain the luminosity at 4392 A rest-frame of each 
source. The adopted average disc inclination throughout this work 
is cos(i) = 0.8 because we select sources that show low obscura- 
tion, and from the AGN standard paradigm are expected to be seen 
face-on or with low inclination. This is also the value adopted in 
DLll, which allows an easier comparison of results. The black hole 
masses are gathered from different sou rces in the following order 
of preference: rever beration mapping jDennev et al.ll2006l 



Peterson et al.l|2004|). Ha vir ial mass estimates ( Winter et al 



201C 



201C 



Ho et al.l l2008': 'Grupe et al."2010l: IParisi et alj|2009h and Jf-band 



luminosity ( Vasudevan et al. 200^. 

The results are shown in Fig.[T2]in blue and in Table[T] We also 
plot the dependence of mass accretion rate with X-ray luminosity 
in Fig. [13] The results show a trend of increasing mass accretion 
rate with X-ray luminosity, which is expected if we assume that 
higher X-ray luminosities are associated with higher optical lumi- 
nosities. We stress that for these sources, the X-ray and optical/UV 
are simultaneous, and that the effect of variability is negligible in 
this case. 
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AGN 


z 


log Mbu/Mq 


log (vL^) 


log L2_iokcV 


log M 


Mass source 


Notes 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


IRXS J045205. 0+493248 


0.029 


8.04+!!'O; 

— 0.01 


44.08 


43.13 


-1.18 


V09 




2MASX J21 140128+8204483 


0.084 


8.68+!!'°i 


44.43 


44.37 


-1.30 


V09 




3C 120 


0.03301 


7.74+!;'2 
' ' —0.2 


44.17 


43.84 


-0.75 


P04 




3C 390.3 


0.0561 


8.46+°'?^ 

^ ^— 0.1 


44.34 


43.39 


-1.21 


P04 




Ark 120 


0.032713 


o .0+0.05 
'J--Lu_o.06 


44.42 


43.85 


-0.80 


P04 




ESO 548-G081 


0.01448 


7.74+O'OJ 

' ' —0.01 


43.04 


42.55 


-2.43 


V09 




IRAS 05589+2828 


0.033 


o r,r, + 0.01 

o-'^^ — O.Ol 


44.07 


43.39 


-1.38 


W09 




IRAS 09149-6206 


0.0573 


7.84 


45.00 


44.28 


0.39 


P09 




MCG +04-22-042 


0.032349 


709+0.01 

—0.01 


43.72 


43.46 


-1.67 


W09 




Mrk 1018 


0.042436 


9c:+0.02 

0»^0 t\ f-\-i 

—0.01 


43.90 


43.74 


-1.66 


W09 




Mrk 279 


0.030451 


7.54+«-l 

— 0.1 


43.63 


42.65 


-1.35 


P04 




Mrk 352 


0.014864 




42.70 


42.02 


-2.45 


H08 




Mrk 509 


0.034397 


8.16+ni 

— 0.04 


44.27 


43.60 


-1.02 


P04 




Mrk 841 


0.036422 


8.05+n? 

— 0.0.2 


43.88 


43.85 


-1.48 


W09 




NGC 4593 


0.009 


6.99+°f 


42.47 


43.67 


-2.55 


D06 




NGC 5548 


0.017175 


7 64+°-^^ 
'■"^-0.09 


43.10 


42.34 


-2.24 


DIO 




NGC 7469 


0.016317 




43.51 


42.56 


-1.08 


P04 




NGC 985 


0.043143 


8.05tH^ 


43.97 


43.77 


-1.35 


GIO 




SBS 1136+594 


0.0601 


8.00t«:S 


43.95 


43.71 


-1.33 


W09 




SBS 1301+540 


0.0299 


7 r,t; + 0.01 
'•^^-0.01 


43.51 


42.92 


-1.24 


V09 




UGC 06728 


0.006518 


6 71+0 03 
"•'^-0.03 


42.43 


43.99 


-2.32 


W09 




WKK 1263 


0.02443 


'•°'-o.oi 


43.37 


42.92 


-1.88 


V09 





SDSS J002332.34-01 1444.1 


0.483 


9.24 ± 0.03 


45.38 


44.48 


-0.42 


SIO 


SDSS J005709.94+144610.1 


0.172 


9.38 ± 0.02 


44.84 


44.37 


-1.38 


SIO 


SDSS J092809.43+383000.5 


0.498 


9.19 ±0.02 


45.55 


44.09 


-0.13 


SIO 


SDSS J093944.53+392402.8 


0.455 


8.66 ± 0.06 


45.30 


43.49 


0.03 


SIO 


SDSS Jl 11706.39+441333.3 


0.144 


8.77 ± 0.02 


44.54 


43.91 


-1.22 


SIO 


SDSS J142424.21+595300.4 


0.135 


8.69 ± 0.05 


44.38 


43.73 


-1.39 


SIO 


SDSS J142455.53+421407.6 


0.316 


8.47 ± 0.09 


44.93 


44.60 


-0.34 


SIO 


SDSS J163051.74+471118.9 


0.270 


8.43 ± 0.07 


44.86 


44.00 


-0.40 


SIO 



Table 1. Results from the optical/UV and X-ray data on the Swift/BAI sample of local AGN and the bright SDSS quasars. (1) Redshift. (2) Log of central black 
hole mass from different available sources (see column 6), with errors. (3) Optical luminosity at A = 4392 A in erg s~^. (4) Log of 2-10 keV luminosity in 
erg s~^ from fit to 0.3-10 keV regime for the Swift/BAI sample. For the SDSS quasars this value is only an approximation since we assume a powerlaw with 
index F = 1. 9 and no intrinsic absorptio n. (5) Mass accretion rate in Mp,/ yi from Eg. [3] (6) So urce for black h ole mass value: D06 and DIO) Reve rberation 
mapp ing from lDenney et alj|200( ,|2010| P04) Reverber ation mapping fro m ipeterson et alj|2004 W09) Hn fromlwinter et alj|2010l ; H0 8) Hn fromlHo et alj 
I2OO8I : GIO) Ha from lGrupe et alj20l'or P09) Hp from lParisi et alj2009l : V09) iiT-band luminosity from lVasudevan et alj2009l and S10) rShen et alj201ll . (7) 
Symbol * indicates sources with more uncertain values - see text for explanation. 



4.2 Infrared observations 

One option to remove the need of assuming a full SED shape on 
the bolometric luminosity calculation, is to measure the emitted 
reprocessed flux in the IR. With this approach, assuming that the 
disc emission is re-emitted by the obscuring heated material, it is 
possible to integrate the observed IR and the X-ray emission and 
obtain an estimate of the bolometric luminosity. 

Some of the sources in the BAT sample for which we deter- 
mined M, have been observed by the Infrared Ast ronomical Satel- 
hte ( IRAS) and analysed in a thorough way by Vasudevan et al.l 
( I2OIQ ) following the method of IPozzi et ai] ( I2OO7I) . They used the 



X-ray measured absorbing column density to choose an IR emis- 
sion template and corrected their values for the disc covering fac- 
tor, which determines what fraction of the emission will be repro- 
cessed. The host galaxy IR contribution was also corrected using 
two methods which give similar results: the correlat ion found be- 
tween the nuclear 12.3 jxm and 2-10 keV luminosity dOandhi et al.l 
l2009i) . and host galaxy IR SED templates. In the following calcu- 
lations, we adopt the values of integrated I/i-ioooji-m using the 



latter correction method, and calculate the bolometric luminos- 
ity using: Lboi = L(i-iooo) |xm + ix(o.5-ioo)kov (but use the 
iGandhiet al]|2009l method for calculations with high-resolution 
data below). There are more sophisticated methods o f estimating 
the IR reprocessed bolometric output. For example, IPozzi et aU 
( l2010i) use the torus model from iFritz et all JioO^ on higher red- 
shift quasars, but rec ent clumpy torus models such as the one by 
iNenkova et alj ( |2008|) are also able to explain detailed features in 
mid-IR sp ectra. We neverthe less use the results from the simple 
method by'Pozzi et al. (2007) employed in lVasudevan et al] bOlOl) 
here, since they are available for the sample of interest. The accre- 
tion efficiency is then obtained using these bolometric luminosities 
and the M values calculated in the previous section. Our results 
are shown in Table[3]and in Fig.[T4] Some of the sources lie above 
the maximally rotating efficiency limit for a black hole, which lead 
us to a discussion on possible errors associated with this method. 
Since these are systematic errors and not statistical, the first obvious 
question is: are the black hole masses overestimated? This could be 
true for some of the abnormally high efficiency sources, which have 
their masses determined from the A'-band bulge luminosity / black 
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Figure 12. Mass accretion rate as a function of black hole mass for our two 
samples of AGN and the PG quasar sample of DLl 1 (black stars). Green 
squares represent higher redshift bright quasars, and blue circles the low 
redshift SwiftfBPCT AGN sample. Open blue circles indicate sources with 
more uncertain measurements (the SED shows a sfrange shape that could 
be due to uncertain black hole mass or mild presence of dust). Some data 
points have small error bars which are not visible at the scale of the plot. 
The data for this tigure can be found in Table[T] 



Figure 13. Mass accretion rate as a function of X-ray (2-10 keV) luminosity 
for the Swift/BAI local AGN sample. The optical luminosity used to derive 
M and X-ray luminosity are from simultaneous observations. Filled and 
open circles as in Fig. 1121 



hole mass correlation and from virial methods respectively. These 
methods have associated uncertainties, specially the J-f-band cor- 
relation that depends on the bulge determination (AGN and galaxy 
disc subtraction). 

One of the sources clearly above the limit is NGC 5548. It has 
a relatively accurate black hole mass value determined from rever- 
beration mapping, so its high accretion efficiency is most likely 
caused by an overestimated bolometric luminosity. The archival 
IRAS data used in this analysis were taken in 1983, while the 
UVOT data are from a 2007 observation. We are not probing con- 
temporaneous events, which can be important in the presence of 
flux variability. In fact, this source has been shown t o have long 
term optical and X-ray variability (e. g. 1993 to 1999 - Merkulova 
l2002i and 1996 to 2001 - Utfley et aL| r2003l) . We searched the lit- 
erature for optical measurements closer in time to the IRAS data 
campaign and f ound the specific flux at 439 A for this source mea- 
sured in 1974 l lde Bruvn & Sargen3ll978l) . Simply based on the 
light travelling times for radiation coming from different parts of 
the AGN, we expect the IR emission to lag the optical by a fac- 
tor of months up to years depending on the obscuring gas location 
(parsecs to tens of parsecs). In any case, the 1974 flux is almost a 
factor of 1.9 higher than the one from 2007, which means that if 
we repeat our calculations using the 1974 optical flux and the 1983 
IRAS IR flux, we would obtain an accretion efficiency of roughly 
half the value in our table. This would still be above the limit. It 
is possible that other systematic uncertainties effect our measure- 
ments, such as the constant inclination and accretion disc model 
used. The more complex models described in DLIl give higher 
values for the mass accretion rate, and lower values for the accre- 
tion efficiency than our simple relations. This could contribute to 
the higher values observed here. Going back to the flux variabil- 
ity, NGC 5548 is an example of the variations one can obtain when 
using non-contemporaneous data to determine the accretion effi- 
ciency. 

We also note the very low efficiency IRAS 09149-6206. This 



source was discovered due to its IR emission, b ut has not be e n stud - 
ied in great detail. We used the mass value bv |Parisietal] ( l2009h . 
which is Mbh = 7 x lO^Mp) but fr om the Ti'-band its mass is 
Mbh = 3 X 10^ A^q jVasudevan et al . 2009) which is a very large 
difference, and moves the source to the upper right region of Fig. [14] 
(the two extreme values are plotted as red squares). The mass value 
is most likely somewhere in between the two values mentioned 
above, the lower value adopted gives an extremely high mass ac- 
cretion rate compared to any other source in our sample. A set of 
observations in t he mid-IR closer in time to our optical/UV data 
has been done bv lGandhietai] ( l2009l) . They used the VISIR instru- 
ment on the VLT to determine the luminosity at 12 /im (Li2fjm)- 
The high resolution of these observations allowed to measure the 
flux from the core of a sample of objects with very low contami- 
nation from the host galaxy. We use for comparison, their data on 
three of our sources (Mrk 509, NGC 4593 and NGC 7469). We use 
their luminosity and bolometric correction (Eq. 5 of their paper) to 
determine alternative bolometric luminosities. We plot our results 
of bolometric luminosity as a function of black hole mass in Fig. 1151 
and other properties in Table[2] The asterisks in Fig.[T5]indicate the 
value of luminosity from the VISIR data and are connected by lines 
to the correspondent values using the IRAS flux. Red circles repre- 
sent the same IRAS sample as in Fig. [14] This example illustrates 
an alternative to obtain the bolometric luminosity. The difference 
in values can be due to several factors, including the simple bolo- 
metric correction used, but most likely has a contribution from flux 
variability. All the AGN in the Swift/BKI sample may be affected 
by variability. The study of individual accretion efficiencies could 
be extended, providing the other uncertainties are accounted for, by 
selecting a sample with contemporaneous opticalAJV, X-ray and IR 
data, where the IR lags the optical and UV emission by a few years. 
There are new surveys for which data is starting to become avail- 
able, such as AKARI and WISE, that target the IR bands, and could 
be used to reconstruct contemporaneous SEDs for a study such as 
the one described above. 
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log Lboi 


e 


-^Edd 




Notes 




(1) 


(2) 


(3) 


(4) 


(5) 


Mrk509 


45.39 


0.46 


1.4e-l 


7.5e-2 




NGC 4593 


44.22 


1.04 


1.4e-l 


1.7e-2 




NGC 7469 


45.08 


0.26 


7.8e-l 


2.7e-2 





Table 2. Results for sources with VISIR IR data fro mlGandhi et alj2009t (1) Log of observed bolometric luminosity in erg s ^ calculated using the X-ray and 
L\2^m correlation and the bolometric correction in iGandhi et alj2009l . (2) Accretion efficiency calculated from Eq.|4] (3) Eddington ratio (A = ibol/iEdd)- 
(4) Bolometric correction Lopt/^bol- (5) Symbol * indicates sources with more uncertain values - see text for explanation. 



4.3 SDSS quasar sample 

As mentioned in the beginning of Section]?] for high mass AGN at 
the right redshift, it is possible to constrain the accretion disc emis- 
sion peak in the optical/UV wavelengths. In this section, we anal- 
yse a sample of 10 bright quasars with SDSS and GALEX data at 
higher redshifts z ~ (0.13 — 0.5). They have flux measurements in 
the u, g, r, i and z SDSS bands and non-simultaneous in the Near- 
UV and Far-UV bands of GALEX, which are corrected for galac- 
tic reddening. The host galaxy contribution was not removed for 
these sources since they are bright quasars and are expected to out- 
shine the host galaxy. There is also a single X-ray measurement, 
although this is not simultaneous with the UV and optical data. A 
detailed analysis of their SED and spectral indices will be done in 
Vasudevan et al. (in prep). Out of the initial sample, two sources 
were excluded due to signs of dust effects on the SED, which leave 
us with a sample of 8 quasars. 

We determine L^, at the rest-frame wavelength of 4392 A, by 
modelling the SED in that region, based on the seven bands for 
which there is data information. The black hole masses used are 
virial mass estimates from iShen et al.l ilOl lb (column 139 in their 
Table 1). Since the redshifts are higher, we do not extrapolate the 
flux between observed bands but assume a simple model to inter- 
polate the data (see below). This should be spin independent as we 
already discussed in Fig.[T] The quasars were selected to have their 
emission peak at the observed bands mentioned above and we de- 
termine the mass accretion rate for 8 sources but narrow down the 
sample for the accretion efficiency measurements. We require to 
have not only the emission peak at the right energies but also a well 
constrained turnover by the data, which is only observed for 3 of 
the sources with similar shapes to Fig.[TO](see list in Table|4j. For 
the Swift/BKI sources, the black hole mass and low redshift set the 
disc emission peak at higher energies than the ones probed by the 
UVOT bands. In this case, there is a wide range of accretion ef- 
ficiencies that are consistent with the measured M. For the SDSS 
sources, due to their higher mass values, the data points are dis- 
tributed around the disc emission peak and we use this to constrain 
the bolometric luminosity. To obtain an approximate SED shape 
based on the data, we use XSPEC to fit a model to our observa- 
tions. Since we only have one X-ray data point, we cannot con- 
strain the amount of obscuration (hydrogen column density), and 
the slope of the X-ray powerlaw. To get an approximate shape, we 
fix the powerlaw index to a typical value of 1.9, and consider that 
our source is only affected by galactic absorption. This means that 
the bolometric luminosity obtained is most likely a lower limit to 
the real value, due to an underestimated X-ray luminosity. This ef- 
fect is nevertheless small, because we expect most of the emission 
to be in the optical/UV region of the SED. For one of our sources, 
SDSS J092809.43+383000.5, the X-ray flux is an upper limit to the 
real value. We nevertheless use this limit in our calculations due to 
the small effect of the X-ray luminosity on the bolometric luminos- 



ity and the uncertainties in the photon index mentioned above. The 
optical-UV region is modelled by a multi-temperature blackbody 
accretion disc around a Kerr bl ack hole, to a llow blackbody shape 
variations with spin (KERRBB, Li et al.ll200 5). Our final model is 
then [KERRBB + WABS(BKNPOWERLAW)], with a redshift factor 
added manually, and intrinsic reddening taken as zero. The intrin- 
sic SED is then integrated from (10"^ — 100) keV to obtain the 
bolometric luminosity. An example of the best-fit model is shown 
in Fig. [To] We emphasize that our main objective is to constrain the 
shape of the SED with the data points to be able to obtain Lboi, 
and do not use the model output for the spin parameter. One of 
the important factors when modellin g the accretion d isc emission 
is the spectral hardening factor /coi dShimura & Taka hara 1995]), 
which gives the ratio between the colour temperature of the black- 
body used to model the emission and the disc effective temperature. 
A change in this parameter will not only shift the modelle d emis- 
sion peak but change the spectral shape JRoss et alj[l992l) . caus- 
ing degeneracies with the spin determination. Here we assumed a 
constant typical value of /coi = 1.7, but note that previous work 
has shown that fcoi c ould also vary with other physical parameters 
jMerloni et a"i]|2000l) . In the model used, a change in /coi corre- 
sponds to a simple energy shift. This adds uncertainties since, as 
mentioned above, we should be including a spectral shape change 
as well. The value we obtain is just an approximation, valid for our 
assumptions, but we wanted to stress the importance of considering 
more complex accretion disc models that can tackle this issue (ex- 
amples of some of the developments in this field are given in DLl 1 
and references therein). The results are shown in Fig.[T6](green cir- 
cles) and in Table|4] The quasars probe the highest black hole mass 
range, but show a spread in accretion efficiency from the slowly 
rotating to highly rotating black holes. They could also be affected 
by the accretion disc model used to determine the mass accretion 
rate as mentioned in the previous section. 

With this method, we can only measure the accretion effi- 
ciency for sources of relatively high masses, since our data has to 
sample the peak of emission to be able to constrain the bolometric 
luminosity. 



4.4 Constraining the SED observationally 

In the previous sections we referred to different methods of deter- 
mining the bolometric luminosity. Ideally, an accurate luminosity 
should be the result of a fully observationally probed SED. The ob- 
served SED includes the obscuration effects and dust reddening, 
together with a full wavelength coverage, one would also have to 
model the obscuration signatures imprinted on an intrinsic emitted 
SED. 

From theoretical arguments, we would expect to have the en- 
ergy peak of the disc emission depending on the black hole mass, 
with lower mass black holes peaking at higher energies. For the 
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Table 3. Results for sources with IRAS IR data available. (1) Log of observed bolometric luminosity in ergs ^ calculated from i^iR(l — lOOOi^m) + Lx 
(0.5 - 100 keV) jVasudevan etai]|2009l) . (2) Accretion efficiency calculated from Eq. |4] (3) Eddington ratio (A = ^bol/^Edd)' (4) Bolometric correction 
^opt/^bol- (5) Symbol * indicates sources with more uncertain values - see text for explanation. 
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Figure 14. Accretion efficiency as a function of black hole mass for 
SwiftfBPd sources with IRAS IR measurements. Open circles represent 
sources with uncertain mass accretion rate measurements (open circles in 
Fig.[l2]as well). Squares indicate the two possible values for IRAS 09149- 
6206. See discussion in text about extreme sources. 



Figure 15. Bolometric luminosity as a function of mass for sources with 
IRAS measurements and for sources with VISIR data from iGandhi et al] 
I2OO9I Circles represent results for the IRAS sample in Fig. 1141 and asterisks 
the three sources for which there are high-resolution VISIR data. 



highest mass black holes (Mbh ~ 10 A^q) such as our SDSS 
quasars, the optical and UV data can be used to constrain the peak 
turnover, but the Wien type fall-over slope is not possible to deter- 
mine because it is expected to be at the energies where radiation is 
absorbed by obscuring gas (intrinsic and/or the Milky-Way line of 
sight material). Without these data, it is difficult to verify the devia- 
tions from a black-body type decreasing slope, which are expected 
from more complex models of accretions discs (see for example 
Fig. 2 of DLl 1). When we go to lower masses, the opticalAJV data 
probes the Rayleigh-Jeans type part of the disc emission. This al- 
lows to determine the mass accretion rate but does not break the 
degeneracies in terms of black hole spin (such as the case for our 
SwiftfBK^ sources). 

One other option, if we probe even smaller black hole masses, 
would be to look at the low X-ray energies {E <2 keV), to s ample 
the Wien tail of the accretion disc emission. I Yuan etal.l ( [2oToh anal- 
ysed the data on RX J1633-I-4718, an AGN with black hole mass of 



the order (2 - 4) x 10® M©. In this work, they claim that the X- 
ray data are probing the higher energy wing of the accretion disc 
emission. Although there would always be regions more difficult to 
probe observationally, this approach could be used to constrain the 
highest energy range for low mass black hole disc emission. 



5 DISCUSSION 

We have measured the mass accretion rates for 22 local Swi/f/BAT 
AGN and 8 SDSS quasars. The results in Fig. [12] show that with 
these two samples, we are, as expected, probing different black hole 
mass ranges. The M distribution shows a large spread for any of the 
samples, but the local, less luminous AGN tend to extend to lower 
values of mass accretion rate than the SDSS quasars and the PG 
quasar sample of DL 1 1 . The distribution of properties for our sam- 
ples is shown in Fig.[T7] where blue circles represent the 5wi/j/BAT 
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Figure 16. Summary plot of our results. Accretion efficiency for all the 
AGN in this work compared with the distribution found by DLll (black 
stars), as a function of their black hole mass. The red circles represent the 
Swift/BfiH sources with bolometric luminosity determined from IR mea- 
surements, and the green squares are the SDSS sources with bolometric 
luminosity calculated from a modelled accretion disc turnover. The hori- 
zontal dotted lines are from top to bottom, the maximum spin, non-rotating 
and maximum counter-rotating spin for a black hole. 

sources and green squares the SDSS quasars. For comparison we 
over-plot the DLl 1 distribution in black. 

To determine the accretion efficiencies, we selected a 
sub-sample of these 30 AGN for which we could obtain better 
estimates of the bolometric luminosity. Two different methods of 
determining the bolometric luminosity were used: constraining the 
turnover of the quasi-blackbody accretion disc and considering the 
reprocessed emission in the IR, using different sets of available 
archival data. The final results are shown in Fig. [161 The accretion 
efficiencies show a spread betwe en the minimum and max imum 
spin case (0.057 < e < 0.31, iNovikov & Thornelll973 l). with 
the exception of 4 sources which were discussed in Section [431 
They do not show a trend with black hole mass, which is what 
is expected from the discussion in Section [3] Since we are using 
different samples with different selection effects which extend 
to lower I/opt, we expect to see a different distribution in the 
accretion efficiency vs black hole mass plane. The SDSS sources 
have lower accretion efficiencies, and from Fig. [5] could appear to 
be above their Eddington ratio. In fact, from Fig. (Ts] where we 
plot the distribution of bolometric corrections, we can see that 
they are in the group of highest corrections. This means that if we 
simulate their distribution as in SectionjS] the constant Lopt limits 
will correspond to a lower Lboi and a lower Eddington ratio as 
well, as can from the Eddington ratio values in Table |4] In Fig. 1 191 
we plot the distribution of accretion efficiencies as a function of 
their Eddington ratio (A — I/boi/iEdd) for the sources with IRAS 
measurements and the SDSS quasars. We see that in general they 
tend to have moderate to high Eddington ratios (A > 0.01) and 
high accretion efficiencies. This will locate them in the unobscured 
(low a bsorption) AGN region p redicted from AGN evolution 
models ( Raim undo & Fabianl20 09'). which is what we expect since 
these sources, local AGN and bright quasars, were selected to have 
low obscuration. 

As we have been mentioning throughout the paper there are 



Figure 17. Black hole mass and optical luminosity distribution for all 
sources. DLll PG quasar sample (black stars), green squares for SDSS 
quasars and blue circles for Swift/BfiJ^ AGN. Open circles indicate sources 
with more uncertain Lopt and are indicated with a different symbol in Ta- 
bleH 



major uncertainties in the determination of accretion efficiencies. 
Here we summarise the main factors: 

- Selection effects: Our data coverage is limited, and selecting 
sources from different surveys will introduce different flux and vol- 
ume limitations. Targeting large samples allow us to use the same 
method on a large scale and have all sources affected by approx- 
imately the same uncertainties, but will only provide access to a 
limited parameter space. This can cause apparent correlations be- 
tween black hole accretion parameters. In this work we analysed 
different types of sources, but are still far from being able to draw 
conclusions on the general behaviour of AGN. 

- Dust and host galaxy: In our analysis, and also in DLl 1 previous 
work, we assumed E{B — V) = for all our sources. This intro- 
duces errors which mainly affect the bolometric luminosity mea- 
surement, underestimating its value. Although we have used very 
bright quasars and low absorption SwiftfBAJ' AGN and excluded 
sources which have signs of dust in their SED, it is probable that 
every source will have some dust contamination. One option would 
be to use E{B — V) measured spectroscopically, but since dust and 
spin are degenerate in simple accretion disc models, errors on the 
dust estimate could affect the bolometric luminosity by a signifi- 
cant factor. In the case of the IR bolometric luminosities this factor 
is less important, since we are integrating the emission, but there is 
still the host galaxy contamination. The host galaxy will not only 
contribute to the IR SED part but also when measuring the optical 
luminosity. Ignoring the host galaxy emission will cause an overes- 
timated Lopt and M and an underestimated accretion efficiency 

- Variability: We addressed this issue in Section 1421 when we de- 
termined the bolometric luminosity by adding the intrinsic X-ray 
emission and the reprocessed optical/UV emission in the IR. Our 
mass accretion rates were determined using recent observations in 
the optical/UV, but the IR data were gathered ~20 years before. For 
a variable source, such as some of the ones we analysed and dis- 
cussed, the flux variation in time can be very large, which means 
that we are measuring the mass accretion rate and the bolometric 
luminosity at different AGN states. One approach to remove this 
effect would be to get contemporaneous data on the IR and opti- 
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Figure 18. Bolometiic coiTectioiis for optical luminosity at 4392 A. Red 
circles represent local AGN sample detected with 5wi/?/BAT and green 
squares SDSS quasars. Sources show less t han a factor of ten spread in 
bolometric corrections, as one would expect iMarconi et alj2004h . 



calAJV, as we exemplified with the VISIR data. The IR is expected 
to lag the optical emission by a few months up to years, which 
means that the observations would not have to be simultaneous but 
close in time. 

- Black hole mass: This is always an uncertainty in AGN studies, in 
our case it will contribute to every type of measurement we do. In 
our sample there are some black holes with reverberation mapping 
masses, which are the most accurate, but for some AGN we only 
have values from correlation with K-has\A host galaxy luminosity, 
which has larger uncertainties. 

Two other assumptions in this work are the accretion disc in- 
clination (applied in every mass accretion rate determination) and 
the assumed accretion disc model (applied in every mass accre- 
tion rate determination and accretion efficiency for SDSS quasars). 
The inclination has a factor of (cosi)"''''^ in the mass accretion 
rate. We used cos(i) = 0.8 to compare our results with DLll but 
also because in the standard AGN paradigm, sources with low ob- 
scuration (such as ours) have most likely a smaller angle between 
the disc normal and our line of sight. A face-on disc establishes a 
lower limit for M, and an upper limit for the accretion efficiency. 
It is probable that AGN show a distribution of different inclinations 
and an average value approach is not valid. The values presented 
in Table [T] can be easily corrected for different inclination angles. 
The disc model we used is a simplification and previous work have 
shown that accretion discs are likely more complex than our sim- 
ple approximation. This method, described in Section |2] is useful 
to establish the main physical relations, and show almost the same 
trends as the best-fit relations found for more complex disc models 
in DLll. 

The methods explored in this paper can be used to target dif- 
ferent types of AGN using a multi-wavelength data coverage. As 
described in Section [44] black holes of different masses and red- 
shifts will have their emission peak at different wavelength ranges, 
making it possible to constrain some characteristics of their spectra 
in limited wavebands. These observations can be used to comple- 
ment the information we already have, providing the sources of 
uncertainties mentioned above are considered. 



Motivated by the interesting work o f' Davis&Lao3j201lh . who de- 
termined the accretion efficiency distribution for a large sample of 
quasars, we investigated the methods and uncertainties in calculat- 
ing the accretion properties of individual AGN. We explored the 
correlation between th e accretion efficiency and black hole mass 
(e oc M~°-^) found bv lDavis & Laoil j201lh . and conclude that it 
is mainly caused by the selection effects in their sample. In terms 
of the uncertainties, we find that the distribution of calculated ef- 
ficiencies will show an artificial spread, due to possible errors in 
the black hole mass and due to assuming a constant inclination for 
sources which intrinsically have a broad distribution of inclination 
angles. We also exemplify how samples at different luminosities lie 
in the efficiency black hole mass plane, and show that an Eddington 
limited sample (such as A < 1) adds a mass dependent boundary. 

To test the method of determining the mass accretion rate and 
accretion efficiency, we select a low redshift sample of AGN with 
available multi-wavelength data and a sample of SDSS quasars 
with higher mass and luminosity. Our final results show that, with 
simplifying assumptions, it is possible to calculate the accretion ef- 
ficiencies for individual AGN, although the values are largely sen- 
sitive to the unknown parameters and uncertainties. From the cal- 
culations, some AGN give accretion efficiency values above that 
for the maximally spinning case, which is not expected. The effi- 
ciencies are affected by uncertainties in the bolometric luminosity 
calculation, variability, accretion disc model, the black hole mass, 
inclination, dust and host galaxy contamination. The systematic er- 
rors dominate over the statistical errors, which makes it difficult to 
understand the general behaviour of AGN, and to determine corre- 
lations without further observational data on the sources. The meth- 
ods explored in this paper can be applied to other samples in the 
future, providing we are able to determine the parameters more ac- 
curately and constrain the errors involved. We discuss how such a 
study can be improved by observing sources with different prop- 
erties, in particular the black hole mass, since their accretion disc 
emission will be expected to peak at different wavelengths. In the 
future, sources with good wavelength coverage, and especially with 
contemporaneous observations to account for variability, would be 
a good starting sample for such a study. 

The determination of mass accretion rate and accretion effi- 
ciency for individual AGN is of great importance in understanding 
the growth of black holes and their relation with the host galaxies. 
The accretion efficiency is related to the black hole spin parameter, 
which has been studied theor etically as a method to learn the mech- 
anisms of black hole growth |yolonteri et al.| [ 2005l , iKing & Pringld 
l2006l iBerti & Volonteril l200ar If determined for large samples, 
these parameters would not only shed light on why black holes 
show different physical properties but would also be able to con- 
strain models of AGN populations and evolution. 



7 ACKNOWLEDGEMENTS 

SR thanks W. N. Brandt for useful comments on the paper and for 
providing the data on the SDSS quasars, A. Zoghbi for kindly help- 
ing to update XSPEC models and R. Canning for useful discus- 
sions. SR acknowledges financial support from FCT - Fundajao 
para a Ciencia e a Tecnologia (Portugal). ACF thanks the Royal 
Society for support. 



© 0000 RAS, MNRAS 000, 000-000 



14 S.I. Raimundo, A.C. Fabian, R.V. Vasudevan, P. Gandhi and J. Wu 



AGN 


log Lbol 


e 


-^Edd 


1/fcopt 




(1) 


(2) 


(3) 


(4) 


SDSS J002332.34-01 1444.1 


46.06 


0.17 


5.2e-2 


0.21 


SDSS J005709.94+ 1446 10.1 


45.55 


0.48 


1.2e-2 


0.19 


SDSS J092809.43+38300().5 


46.20 


0.12 


8.3e-2 


0.22 



Table 4. Results from SED fitting for sources with observed turnover (1) Log of observed bolometric (0.0001-100 keV) luminosity in erg s~^ measured 
by integrating the SED wfhich best fits the UV/optical region. (2) Accretion efficiency calculated from Eq. |4] (3) Eddington ratio (A = Lbol/^'Edd)- (4) 
Bolometric correction iopt/^bol- 
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Figure 19. Accretion efficiency as a function of Eddington ratio for local 
and high redshift AGN (colour coded as in Fig. [let. Our AGN he in the 
region of unobscured sources of Fig. 3 in lRaimundo & Fabianll2009l as ex- 
pected due to our selection criteria of low obscuration sources. 
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APPENDIX A: SPECTRAL ENERGY DISTRIBUTIONS 
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Figure Al. Observed-frame spectral energy distributions for sources in Table[T] Lines are plotted to indicate possible shape of SED (dashed lines for observed 
and solid lines for intrinsic before absorption) assuming E(B — y) = 0. These are not used in the calculations except for SDSS J()02332. 34-011444.1, 
SDSS J005709.94-H144610.1 and SDSS J092809.43-^383000.5. Mass accretion rate is determined direcdy from data points and bolometric luminosity from IR 
measurements for the majority of sources (see description in Section|4]for details). For most of the sources plotted the optical/UV turnover is not constrained 
by the data. We assume a value of spin for plotting purposes. 
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